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ABSTRACT 

We demonstrate a new procedure to derive accurate and precise surface gravities from high 
resolution spectra without the use of external constraints. Our analysis utilizes Spectroscopy 
Made Easy (SME) with robust spectral line constraints and uses an iterative process to miti¬ 
gate degeneracies in the fitting process. We adopt an updated radiative transfer code, a new 
treatment for neutral perturber broadening, a line list with multiple gravity constraints and 
separate fitting for global stellar properties and abundance determinations. To investigate the 
sources of temperature dependent trends in determining logg noted in previous studies, we 
obtained Keck HIRES spectra of 42 Kepler asteroseismic stars. In comparison to asteroseismi- 
cally determined log 5 our spectroscopic analysis has a constant offset of 0.01 dex with a root 
mean square (RMS) scatter of 0.05 dex. We also analyzed 30 spectra which had published 
surface gravities determined using the a/i?* technique from planetary transits and found a 
constant offset of 0.06 dex and RMS scatter of 0.07 dex. The two samples covered effective 
temperatures between 5000K and 6700K with logg between 3.7 and 4.6. 

Subject headings: stars: fundamental parameters; stars: solar-type; techniques: spectroscopic; 
asteroseismology; methods: data analysis 


1. INTRODUCTION 

For many planets, uncertainty in the stellar ra¬ 
dius is the largest contributor to uncertainty in 
planet radius and density, limiting attempts to un¬ 
derstand planet formation theory. Ten years ago, 
the high quality Spectral Properties of Cool Stars 
(SPOCS) catalog (Valenti & Fischer 2005, hereafter 
VF05) provided a uniform analysis which enabled 
new insights into planet formation Fischer & Valenti 
(2005). However, the authors cautioned that offsets 
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remained with other catalogs and that a temperature 
dependent bias existed between the spectroscopic and 
model isochrone gravities for some stars. 

More recent analyses have used non spectroscopic 
(’external’) constraints on either Teff or logt; to re¬ 
move this trend. Torres et al. (2012) showed that the 
Spectroscopy Made Easy (SME) (Valenti & Piskunov 
1996) analysis used in VF05 displayed degeneracies 
between and \ogg and that the surface gravity 
could be off by up to 0.5 dex in regions of the HR 
diagram where the Mg b triplet is a poor gravity 
constraint. An offset of 0.5 dex in logt; resulted in 
an a corresponding offset of 400K in Teff- One re¬ 
cently developed external constraint for stars with 
transiting planets, used in Torres et al. (2012), is the 
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stellar density derived from the a/i?* ratio (Sozzetti 
et al. 2007). Although transiting systems are fan¬ 
tastic laboratories, only 1% of planetary systems will 
be favorably aligned for transits. Additionally, this 
method can suffer from inaccuracies due to contam¬ 
ination from the light of nearby companions (Seager 
& Mallen-Ornelas 2003) as well as high impact pa¬ 
rameter or large eccentricity of the transiting planet 
(Huber et al. 2013). Angular diameter measurements 
of stars are a gold standard for obtaining effective 
temperatures when the angular resolution of the in¬ 
terferometer can directly measure the angular stellar 
diameter and accurate distances are known. This is 
an important means of validation but is not possible 
for more distant stars. 

Valenti et al. (2009) externally constrained SME 
by using parallax to determine the bolometric lumi¬ 
nosity of the star then combined that with derived 
spectral parameters to interpolate in a grid of stel¬ 
lar models and obtain a constraint on gravity. Gaia 
will eventually give us precise distances for millions 
of stars; however, many of the Kepler planet hosts do 
not have well-measured distances. 

To avoid using external constraints, we have 
searched for the source of the inaccurate gravity de¬ 
terminations and correlated errors in temperature 
and metallicity. Most of the gravity information in 
the line list of VF05 is contained in the damping 
wings of the Mg I b triplet lines, which are sensitive to 
pressure changes in main sequence stars cooler than 
about 6200 K. The visual spectrum is dominated by 
lines of Fe I which is the minority species in the Sun 
and hotter stars. With few Fe II lines in the line list 
the solution is less sensitive to the ionization equilib¬ 
rium. Our initial investigation showed that although 
there were minor improvements which could be made 
by improving the treatment of pressure broadening 
there was still a large degeneracy between temper¬ 
ature, gravity, and metallicity which could lead to 
inaccurate results. 

The line list of VF05 contains ~ 1000 atomic lines 
spanning 170 A and includes a handful of the stronger 
molecular lines and 13 Fe II lines. The upgrade of the 
HIRES CCD from one to three detectors extended the 
wavelength range and we take advantage of this in 
our new analysis. Our line list now contains roughly 
7500 lines covering more than 350 A. The addition of 
new temperature and pressure sensitive lines, includ¬ 
ing 290 Fe II lines, has helped to break the degen¬ 
eracies between temperature, gravity, and metallic¬ 
ity. In addition to the large line list, we also use the 
prescription for broadening by neutral hydrogen from 


Barklem & O’Mara (1998) when available in VALD-3. 
This more accurate treatment of line wing broaden¬ 
ing provides a better fit to the wings of strong lines. 
Through two independent comparisons, we show that 
the procedure presented in this paper makes system¬ 
atic trends in derived surface gravities comparable to 
or smaller than random errors. 

2. OBSERVATIONS AND REDUCTIONS 

Our purpose in this work was to improve spectral 
synthesis modeling for accurate determination of ef¬ 
fective temperature and surface gravity. Kepler as- 
teroseismic observations (Borucki et al. 2010) allow 
precise determinations of stellar mass and radius giv¬ 
ing both accurate and precise surface gravities (Gai 
et al. 2011; Basu et al. 2010). We obtained Keck 
HIRES spectra for a sample of 42 stars covering a 
range of temperatures, surface gravities, and activity 
levels to explore how changes in the spectral anal¬ 
ysis affected our ability to recover the asteroseismic 
surface gravity. 

2 .1. Spectra 

We obtained 43 spectra of 42 Kepler asteroseismic 
targets using the Keck HIRES spectrograph in the 
red conhguration at a resolution of i? ~ 70,000. The 
signal-to-noise ratio (SNR) was typically > 200 per 
pixel column in the region around the Mg I b triplet 
for all but 4 fainter targets that had S/N of ~ 120. 
All spectra were reduced using the standard pipeline 
of the California Planet Search team. 

3. ASTEROSEISMIC SURFACE GRAVITIES 

Stellar oscillation frequencies were obtained using 
short cadence observations from the primary Kepler 
mission (Borucki et al. 2010). The short cadence and 
precision photometry of the Kepler telescope allows 
Aiy and Vmax to be determined extremely accurately 
for most stars (Chaplin et al. 2014). 

To close approximation the large frequency separa¬ 
tion, Al/, scales as where p is the mean stellar 
density and the frequency of maximum power, i/max-, 
scales as (Christensen-Dalsgaard et al. 2010). 

Combining these parameters with an externally de¬ 
termined effective temperature, we can use the scal¬ 
ing relations to find the surface gravity. The uncer¬ 
tainties in logg can be reduced by using a grid of 
stellar evolutionary models with estimated frequen¬ 
cies instead of simply using the scaling relations. We 
used the grid-based Yale-Birmingham pipeline (Gai 
et al. 2011; Basu et al. 2010) which combines Av and 
i^max derived from Kepler lightcurves with the spec¬ 
troscopically determined TeS and [Fe/H]. 
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3.1. Iterative Fitting 

The asteroseismic surface gravities of main se¬ 
quence stars depend only weakly on the effective tem¬ 
perature and metallicity of the star (Gai et al. 2011) 
with their importance increasing slightly as the star 
evolves. We obtained temperature and metallicity 
from our spectral analysis, then used those values 
as initial inputs in the asteroseismic analysis. We 
then iterated once, fixing the gravity in our spectral 
analysis to the value returned from asteroseismology 
and used the newly derived Tes and [Fe/H] values 
in the asteroseismic analysis. The iteration resulted 
in changes of less than 0.05 dex in the final gravities 
from those initially determined from the asteroseis¬ 
mic analysis. 

4. SPECTROSCOPIC ANALYSIS TECHNIQUE 

SME combines a stellar atmosphere grid, an atomic 
and molecular line list, and a radiative transfer code 
to create model spectra based on specified physi¬ 
cal parameters such as effective temperature, sur¬ 
face gravity, metallicity, and rotation. Additionally, 
SME can fit an observed spectrum using Levenberg- 
Marquardt least-squares fitting with any number of 
free global stellar parameters combined with zero or 
more free elemental abundances. 

VF05 analyzed nearly 2000 spectra using 1000 
lines over ~ 170 A to produce the SPOCS catalog. 
Initial line parameters were obtained from the Vienna 
Atomic Line Database (VALD) (Kupka et al. 2011). 
Then VF05 tuned line position as well as log{gf) 
and van der Waals broadening coefficients to fit a 
solar atlas. In addition, 78 strong molecular MgH 
and C2 lines from (Kurucz 1993) were also included. 
Their analysis solves simultaneously for the global pa¬ 
rameters surface gravity (log g), effective temperature 
(Teff), metallicity ([M/H]), projected rotational ve¬ 
locity (?;sini), and radial velocity {vmd)- In addition 
to the global parameters, VF05 solved for elemental 
abundances for Na, Si, Ti, Fe, and Ni. Microtur¬ 
bulence was fixed at 0.85 km/s and they derived an 
empirical relation for macroturbulence as a function 
of effective temperature. The resulting parameters 
had good relative precision, though for some stars the 
spectroscopically determined gravities were inconsis¬ 
tent with isochrone gravities, especially for warmer 
stars. This later motivated the use of external con¬ 
straints on surface gravity in a more recent analysis 
to reduce parameter degeneracy (Valenti et al. 2009). 

4.1. Baseline Analysis with New Code 


As a baseline, we began by analyzing the spectra 
using the line list of VF05 with the same version of 
SME they used. We then updated the line list over 
the same wavelength region and updated the version 
of SME to v439^, which uses an updated radiative 
transfer code and updated atmosphere interpolation 
algorithm. These changes reduced the differences in 
spectroscopic gravities relative to our asteroseismic 
reference values by about half; however, trends in 
these differences as a function of both metallicity and 
temperature remained. Between 5000 K and 6500 
K, A log 5 (spectroscopic minus asteroseismic log 5 ) 
spanned 0.5 dex. 

We then included additional spectral intervals and 
decreased the number of simultaneously free parame¬ 
ters in our model to decouple the fitting of global pa¬ 
rameters from individual abundances. These changes 
dramatically improved the accuracy in our derived 
logg values, reducing the RMS scatter in A \ogg to 
0.1 dex and removing the trend with respect to Teff. 
Further improvements were achieved by iterating the 
fitting procedure using the derived abundance pat¬ 
tern and allowing the abundances of alpha elements 
to be free while fitting the global parameters. This 
approach removed the trend in derived gravity with 
respect to [M/H] and reduced the RMS scatter in A 
logg to only 0.05 dex. 

4.2. Expanded Line List 

The line list now includes more than 350 A in 20 
segments between 5160 A and 7800 A and includes 
nearly 7500 atomic and molecular lines (Table 1). 
The lines were obtained from the VALD-3 database 
(Kupka et al. 2011) and then tuned to better match 
a high resolution disk-integrated solar atlas (Wallace 
et al. 2011). The VALD-3 data contains many as- 
trophysically tuned line parameters, but we found 
that about 15% of the VALD-3 lines still needed ad¬ 
justment in one or more of line position, log(g/), 
or neutral perturber broadening coefficients in order 
to match the solar spectrum. Where available, we 
also used VALD-3 values for the temperature depen¬ 
dent broadening prescription of Barklem & O’Mara 
(1998), which gives better fits to the line profiles than 
the traditional van der Waals formulation. 

We added wavelength regions to increase the num¬ 
ber and quality of gravity and temperature dependent 
lines and to allow abundance determinations for ad¬ 
ditional elements. SME calculates from the differ- 

^ Release versions of SME can be downloaded from 
http;//www.stsci.edu/-valent1/sme.html 


4 



Figure 1. We increased the amount of Tgff sensitive regions 
by 51% and logg sensitive regions by 28% by including addi¬ 
tional wavelength segments. The dotted red line are for those 
segments which correspond to the original wavelength range of 
VF05 and the solid blue line corresponds to the newly added 
wavelength coverage. The gray vertical lines denote the bound¬ 
aries of the wavelength segments and the letters at the top 
denote are keys to the wavelength ranges listed in Table 1. 

ence between observed and model spectra. To evalu¬ 
ate the influence that the new segments had on fitting 
for logg and Teff we created a model spectrum with 
solar parameters, two models with logg adjusted by 
±0.2 dex, and two adjusted by ±200K in Teff. We 
differenced pairs of models and used the squared dif¬ 
ference at each unmasked spectrum point as a proxy 
for the information contributed to when fitting 
to spectra. We then plotted the cumulative distri¬ 
bution of these squared differences to examine this 


gravity and temperature information as a function 
of the increasing wavelength point in our unmasked 
spectrum (Figure 1). The information content of 
each segment is also detailed in Table 1. The ex¬ 
panded line list added 28% more gravity information 
and 51% more temperature information. Addition¬ 
ally, although deep gravity sensitive lines tend to get 
stronger with increasing logg and temperature sen¬ 
sitive lines stronger with decreasing Teff, there are 
some lines which display the opposite behavior. As 
noted by Gray (2008), weak lines of ions or atoms of 
an element that is predominantly found in the same 
ionization state (e.g. weak Fe II lines) provide im¬ 
portant gravity information because these lines be¬ 
come stronger with decreasing gravity. The opposing 
line growth of these lines is especially helpful in con¬ 
straining Teff and logg and we more than doubled 
the number of these spectral lines. Since Fe I is in 
the next lower energy state than the dominant Fe II, 
weak lines of Fe I will not be very pressure sensitive 
and so ionization equilibrium provides an additional 
gravity constraint to the Mg I b wings, for lines where 
non-LTE effects do not significantly affect ionization 
equilibrium. 

In preparing the spectra for analysis we use the 
same continuum normalization procedure as VF05 
but also mask out telluric lines as our spectral range 
now includes regions with significant telluric contam¬ 
ination. We perform cross correlation with the solar 
atlas (Wallace et al. 2011) to find an approximate ra¬ 
dial velocity and before shifting the spectrum to ob¬ 
servatory wavelengths, we apply a telluric mask based 
on telluric lines found in the Wallace et al. (2011) 
atlas. This masking had a beneficial effect on con¬ 
tinuum normalization in some segments for spectra 
where telluric lines fell on continuum features near 
the ends of the segment. 

4.3. The New Analysis Procedure 

Initial stellar parameters are all set to solar val¬ 
ues except for the temperature, determined by V-K 
or B-V color relation (Boyajian et al. 2013), and 
gravity which is arbitrarily set to 4.5. In the first 
step (“Step 1”), we allow the global parameters Teff, 
logg, [M/H], and macroturbulence (umac) to be free 
along with individual abundances for Ca, Si, and Ti. 
These alpha elements represent the largest number of 
non iron peak lines in the spectral regions we are an¬ 
alyzing. They also have relatively high abundances 
that can differ greatly from the solar abundance pat¬ 
tern especially for low metallicity stars. In our ini¬ 
tial fits, overall metallicity is strongly correlated with 
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Segment 

^start 

^end 

Decreasing logp 
lines grow lines weaken 

Decreasing Teff 
lines grow lines weaken 

VF05 

a 

5164 

5190 

0.17% 

67.07% 

48.14% 

0.00% 

X 

j 

5190 

5207 

0.63% 

10.20% 

16.17% 

0.00% 


k 

5232 

5262 

1.93% 

5.77% 

10.78% 

0.01% 


b 

6000 

6015 

0.08% 

0.33% 

1.68% 

0.02% 

X 

c 

6015 

6030 

0.02% 

0.78% 

2.20% 

0.00% 

X 

d 

6030 

6050 

0.29% 

0.03% 

0.35% 

0.03% 

X 

e 

6050 

6070 

0.10% 

0.35% 

1.08% 

0.02% 

X 

f 

6100 

6120 

0.21% 

0.65% 

1.88% 

0.01% 

X 

g 

6121 

6140 

0.30% 

2.92% 

3.89% 

0.00% 

X 

h 

6143 

6160 

0.60% 

0.14% 

1.11% 

0.04% 

X 

i 

6160 

6180 

0.04% 

4.09% 

5.42% 

0.00% 

X 

1 

6295 

6305 

0.11% 

0.67% 

0.78% 

0.00% 


m 

6311 

6320 

0.02% 

0.11% 

1.34% 

0.01% 


n 

6579 

6599 

0.16% 

0.09% 

1.13% 

0.04% 


o 

6688 

6702 

0.00% 

0.01% 

0.22% 

0.00% 


p 

6703 

6711 

0.00% 

0.00% 

0.43% 

0.00% 


q 

6711 

6718 

0.01% 

0.09% 

0.63% 

0.00% 


r 

7440 

7470 

0.16% 

0.55% 

1.16% 

0.01% 


s 

7697 

7702 

0.00% 

0.30% 

0.20% 

0.00% 


t 

7769 

7799 

0.64% 

0.36% 

0.53% 

0.68% 



Table 1 

Wavelength ranges of spectral regions used in this analysis. Those with a check in the VF05 column cover the same wavelengths 
as segments from VF05. The letters correspond to those at the tops of the plots in Figure 1. The columns for decreasing logg 
and Teff quantify the total amount of information contributed to the determination for each segment as plotted in the figure. 


[Fe/H] due to the preponderance of iron lines. By let¬ 
ting these alpha elements be independent of overall 
metallicity, we partially account for non-solar abun¬ 
dance patterns in fitting the global parameters. VF05 
showed that microturbulence, important in equiva¬ 
lent width analysis, seems not to play a large role in 
forward modeling and is degenerate with [M/H] so 
we hxed it to our adopted solar value of 0.85 km/s. 
After fitting, we perturb the temperature ±100K in 
the resulting model and re-fit, using these new initial 
parameters to encourage the non-linear least squares 
solver to converge to a nearby minimum. The new 
analysis provides more consistent results with median 
standard deviations between the three models of UK, 
.01 dex, and .005 dex for Tes logg, and [M/H] for the 
stars in our sample. 

In the next step (“Step 2”) we fix the global param¬ 
eters to the weighted average of the three models 
from Step 1. We then allow all of our elemental abun¬ 
dances to be free including the three alpha elements 
which were free in Step 1. We have a total of 15 
free elemental abundances (C, N, O, Na, Mg, Al, Si, 
Ca, Ti, V, Cr, Mn, Fe, Ni, Y), which were selected 
based on our interest in the element, the number of 
lines available in our spectral range, the sensitivity 


of those lines to small changes in abundances, and 
the reliability of recovering the solar abundances in 
our asteroid spectra (see below). A few of those ele¬ 
ments, carbon, nitrogen, oxygen, and magnesium are 
also components of molecular species with lines in 
our spectral range and our inclusion of molecules in 
the spectral synthesis model provide additional con¬ 
straints on these elemental abundances. The resul¬ 
tant model now has the full set of global parameters 
and individual abundances. 

We then iterate, repeating Step 1 and Step 2 using 
the abundance pattern derived from the output of the 
first iteration instead of solar values. The differences 
in the models between the first and second iterations 
are small but serve to erase a subtle trend in A log g 
as a function of metallicity that was seen at the end 
of the first iteration. The standard deviation of the 
differences between the first and second iteration is 
0.056 dex in logg and only 0.01 dex for [M/H]. 

To evaluate the accuracy of our returned param¬ 
eters, we analyzed 20 spectra of reflected sunlight 
from 4 different asteroids taken on 6 different epochs. 
These spectra were obtained in the same manner as 
our other targets. Because we calibrated the atomic 
line data to the high resolution NSO solar spectrum. 
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we expect that we should recover the solar parameters 
when we fit the reflected solar spectra in the asteroid 
observations. We found that our procedure recovers 
the solar parameters for the asteroid spectra (Figure 
2) with weighted mean offsets of only IK in Teff, 
0.02 dex in log 5 , and 0.02 dex in [Fe/H] and RMS 
scatter of 5K, 0.006 dex, and 0.003 dex respectively. 



Figure 2. The wings of the Mg I b triplet contain significant 
gravity information in our spectral region. As can be seen in 
this fit (heavy red line) to a spectrum from the asteroid Vesta 
(thin black line), we obtain a good fit using our new analysis. 
Gray regions are excluded from the fit and residuals (observed 
- model) are shown below. 

4.4. Microturbulence 

Our analysis is one of differential solar measure¬ 
ments made after tuning our line list against a solar 
atlas. In doing so, we adopted a microturbulence 
value of 0.85 km/s for the sun and left that fixed in 
the analysis of other stars. After fitting using our 
new procedure, we analyzed the effects of this deci¬ 
sion to follow VF05 in fixing the microturbulence to 
the value we adopted for the sun with two tests. The 
first fixed Umic using the final parameters from our 
analysis in the empirical formula of Ramirez et al. 
(2013) 

rmic ~ 1.163 

-h7.808 X 10-'‘ X (Teff - 5800) 

-0.494 X (logg-4.30) 

-0.05 X [Fe/H] 


and then ran our analysis as before. For the aster¬ 
oid spectra Umic = 1-07 was 0.2 km/s higher than our 
value of 0.85 and resulted in models on average 30K 
hotter and 0.03 dex higher in logg. For the astero- 
seismic stars, the resultant gravities were discrepant 
from asteroseismic values by up to 0.35 dex and there 
were clear trends in A log 5 with both Tes and [M/H]. 
Temperatures also increased for hotter stars by up to 
200K. Allowing for the possibility that the offset in 
Vmic at solar from our adopted value could be re¬ 
sponsible for the differences, we re-ran this test after 
subtracting a constant 0.22 km/s from the rimic val¬ 
ues. By construction, this then returned the solar 
parameters for the asteroid spectra but the trend in 
A log g with Tes remained with a slightly shallower 
slope (Alogg up to 0.23 dex). 

The second test allowed Vmic to be an additional 
free parameter in the global parameter step with its 
initial value set to 0.85 km/s and our procedure was 
then run as before. There was a systematic offset 
of ^ 0.2km/s lower in the returned Umic with re¬ 
spect to the fixed test values. The asteroid spectra 
returned the solar parameters with only [M/H] in¬ 
creasing slightly by 0.01 dex. However, asteroseismic 
stars again showed clear trends in A log g with both 
Teff and [M/H] though slightly less (A log g up to 0.25 
dex) than the fixed Umic case. 

The first step in our analysis was to empirically 
tune our line data against a solar atlas using fixed 
solar parameters including 0.85 km/s for microtur¬ 
bulence. Using the formula of Ramirez et al. (2013) 
for microturbulence when fitting stellar spectra de¬ 
graded our gravity accuracy because our lines had 
been tuned at this constant, lower value. It is possible 
that self-consistent use of the Ramirez et al. (2013) 
formula to tune the line data and fit stellar spectra 
might improve the accuracy of our abundance deter¬ 
minations without compromising the accuracy of the 
gravity. 

In our models, microturbulence is partially degen¬ 
erate with \ogg, Teff and [M/Hj. Including it as an 
additional free parameter reduces our ability to re¬ 
cover accurate surface gravities. Solving for micro¬ 
turbulence while tuning the atomic line data might 
allow us to solve for microturbulence when fitting 
spectra without compromising the fit for solar-type 
stars. From these tests it is clear that we have to 
leave Umic set at the value we adopted for the sun 
since this is the value used in tuning our lines against 
the solar atlas. 


( 1 ) 


5. RESULTS 
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Figure 3. The unconstrained values of logg determined by 
our new analysis are much closer to the asteroseismic values 
than those returned by our analysis using the VF05 Line list. 
Moreover, there are no trends with temperature and compar¬ 
isons with gravity and metallicity also show no trends. The 
lone outlier in our new analysis comes from a star with total 
rotational broadening (urot) greater than 25 km/s. 

Comparing all of the stellar parameters derived 
from spectral analysis to known stellar parameters 
is only possible for the Sun, where we have indepen¬ 
dent methods for determining age and composition. 
Even then, debate continues surrounding the accu¬ 
racy of abundance patterns. However, for individual 
parameters such as surface gravity, we now have two 
methods that provide consistency checks for our anal¬ 
ysis. 

5.1. Comparison to Asteroseismic Surface Gravity 

The first of these is to compare our spectroscopic 
gravity with the surface gravity from the grid based 
asteroseismic method described in § 3 as our fiducial. 
Our new analysis process results in surface gravities 
which match the asteroseismically determined values 
of the stars in our sample with an RMS scatter of 
0.05 dex and with virtually no offset (0.01 dex) in 
the zero-point. Improvements with respect to the 
prior analysis using the VF05 procedure can be seen 
in Figure 3. There are no significant trends in A 
log(; between the asteroseimically and spectroscopi¬ 
cally determined values with respect to Teff, logg, or 
metallicity (Figure 4). Errors in these three param¬ 
eters have been shown to be correlated in previous 
spectral synthesis modeling analyses using the VF05 
line list and an older SME version (Torres et al. 2012). 
The extreme outlier at A logg = -0.2 dex (Figure 4) 
is a star with total rotational broadening (v sin i and 


^^mac) greater than 25 km/s which smooths away most 
of the information needed to determine gravity and 
temperature. 

5.2. Comparison with Published Results 

A second comparison of our spectroscopic analysis 
comes from the gravity constraint of transiting ex¬ 
oplanets. We compared our spectroscopic logg val¬ 
ues with those derived by Torres et al. (2012) using 
the a/i?* density method (Sozzetti et al. 2007) for 
stars with transiting planets. Most of the spectra 
have SNR> 100 per wavelength bin, though a quar¬ 
ter of them have SNR below 100. Torres et al. (2012) 
found systematic trends between their gravities based 
on a/R^, and SME gravities obtained using the VF05 
procedure. In contrast, we find no systematic trends 
between Torres et al. (2012) gravities and our SME 
gravities obtained using more spectral segments and 
a two-stage htting procedure. Our SME methodology 
works better than the VF05 procedure for the Torres 
et al. (2012) stars, which are generally warmer than 
the original VF05 sample. We find a constant offset 
of 0.06 dex in logg, relative to Torres et al. (2012). 
Such a small offset could be a result of errors in our 
spectroscopic analysis or may be the result of small 
errors in the a/R* analysis that arise because of in¬ 
accuracies in the impact parameters or eccentricities 
(Huber et al. 2013). 

6. DISCUSSION 

We believe that the largest improvement in our de¬ 
termination of surface gravity was a result of adding 
new wavelength segments. The expanded line list 
added new gravity and temperature dependent lines 
that helped reduce parameter degeneracy. We also 
included a factor of 10 more iron lines in the mask, 
bringing the total to ^ 900 with almost one third 
of them Fe II lines. All of the iron lines, regardless 
of ionization state or excitation potential, are mod¬ 
eled with the same iron abundance. Because most 
of the iron is in an ionized state at the temperatures 
we are interested in, weak Fe II lines will become 
stronger with decreasing gravity while weak Fe I lines 
are largely gravity insensitive. Having a large num¬ 
ber of lines in both states provides another important 
constraint on the gravity and helps to further limit 
the parameter space. The addition of the expanded 
wavelength range gave us very accurate gravities and 
left us with only a small residual trend with metal¬ 
licity. 

The wings of the Mg I b triplet lines provide a 
wealth of information on the surface gravity of a cool 
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log(g) from This Work Tg„ from This Work [M/H] from This Work 

Figure 4. The spectroscopically determined gravities return the asteroseismically determined gravities with only 0.05 dex RMS 
scatter and show no trends with logg, Teff, or [M/H]. In the first panel there is no trend in derived gravity for subgiants through 
solar type dwarfs. In the second panel, we see no trends in derived gravity for stellar temperatures between early K and late F 
objects. In the final panel we see no trends for stars with metallicities between -0.5 and 0.3 dex. The outlier at -0.2 dex is a star 
with total rotational broadening greater than 25 km/s 





Figure 5. A comparison of stars analyzed by Torres et al. (2012) made use of a/R* for transiting planets. There does not appear 
to be any trend in derived surface gravity with respect to the parameters log 5 -, Tefr or [Fe/Hj. This sample contains several low 
mass stars missing in the asteroseismology sample. Spectra with SNR < 100 per wavelength bin have square markers. 
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dwarf star, but the gravity information is degenerate 
with the Mg abundance. This led us to separate the 
global parameter fitting from the abundance determi¬ 
nations. However, because the particular atmosphere 
used is tied to the metallicity, a scaled solar abun¬ 
dance pattern (dominated by [Fe/H]) will not work 
for all stars. By allowing the most abundant alpha 
elements (minus Mg) to be independent of the overall 
metallicity we obtained a closer match to the correct 
atmosphere on the first iteration. The magnesium 
abundance, however, will be the scaled solar value in 
this step and not necessarily accurate. In the sec¬ 
ond iteration, we have a direct determination of the 
magnesium abundance and so obtain more accurate 
gravity information from the Mg I b wings. In Fig¬ 
ure 6 we see that the changes in the magnesium to 
metallicity ratio between the first step (scaled solar 
value) and the final (independently determined Mg 
abundance) are inversely proportional to the change 
in gravity. 



Figure 6. During the first iteration of our analysis, the Mg 
abundance is just the solar abundance pattern scaled to the 
overall metallicity. After obtaining a close estimate of the 
global parameters we get more detailed abundances tor 15 el¬ 
ements including Mg and this new abundance pattern is used 
in fitting for the final global parameters. Because Mg lines 
constitute the bulk of the surface gravity information for cool 
dwarfs, this iterative process allowed us to remove subtle trends 
in derived log(g) with metallicity in our sample. 

As K and M dwarfs are some of the most promis¬ 
ing candidates for finding Earth-like planets in the 
habitable zone a robust spectroscopic analysis of cool 
dwarfs is an important improvement to the work pre¬ 
sented here. We are currently finishing the analysis 


of all stars observed by the California Planet Search 
since late 2004 using this procedure and will include 
our updated line list in that paper. 

7. CONCLUSIONS 

Accurate stellar parameters can be obtained from 
high resolution spectroscopic analysis when spectral 
segments contain adequate constraints and care is 
taken to decouple degeneracies in the fitting process. 
The results improve upon our past analyses where 
correlated errors led to surface gravity determination 
off by more than 0.3 dex for stars where Mg I b is 
no longer a good gravity constraint. Our technique 
returns gravities consistent with those determined by 
asteroseismic analysis with an RMS scatter of only 
0.05 dex. It is also consistent with gravities deter¬ 
mined using a/i?* constraints from planetary transits 
with a small systematic offset of 0.06 dex. 
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